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Abstract 

We perform a quantization of 4-dimensional Nambu-Goto the- 
ory of open string in light cone gauge, related in Lorentz-invariant 
way with the world sheet. This allows to obtain a quantum the- 
ory without anomalies in Lorentz group. We consider a special 
type of gauge fixing conditions, imposed in oscillator sector of the 
theory, which lead to a relatively simple Hamiltonian mechanics, 
convenient for canonical quantization. We discuss the algebraic 
and geometric properties of this mechanics and determine its mass 
spectrum for the states of spin singlet 5 = 0. 



1 Introduction 

This paper continues the work \ : , further referred as Part I, 
devoted to Dirac's quantization of Nambu-Goto theory of 
open string, formulated in the space-time of dimension d=4. 
The approach consists in the introduction of light cone gauge 
with the gauge axis related to a certain dynamical light-like 
vector in the theory. In this approach Lorentz group trans- 
forms the world sheet of the string together with the gauge 
axis and is not followed by reparametrizations. It is shown 
in Part I, that under these conditions the quantum Lorentz 
group is free of anomalies. The details of the mechanics 
depend on the choice of the dynamical vector to which the 
gauge axis is attached, i.e. on the choice of gauge fixing con- 
ditions. Generally the mechanics appears to be algebraically 
complex. In this paper, Section 2 we consider gauge fix- 
ing conditions, producing a resolvable system of polynomial 
equations. We discuss the structure of the obtained solu- 
tions, including their remaining discrete gauge symmetries. 
In Section 3 we perform quantization of this mechanics and 
determination of its mass spectrum for the states of spin 
singlet S = 0. In appendices we present the formulae for 
the solutions of polynomial system, provide proofs of the 
statements formulated throughout the paper and give the 
details about numerical methods used for determination of 
mass spectrum. 



2 Classical mechanics 

2.1 Lorentz-invariant abelian light cone gauge (lia-lcg) 

We start from the canonical basis of string mechanics con- 
structed in Part I: 



PfiiZp + infinite set of oscillators a k , a k + the top e;, S, 
with Poisson brackets 



{a k ,a* n } = ikS kn , fc,neZ/{0}, 
{S\S 3 } = -e l]k S k , {S\ei} = - 



and equivalent symplectic form 
Q, = dP u A dZ 



da* k A da k + \d&\ A d(S x e»). (2) 

fc#0 



The mechanics is restricted by 4 constraints of the 1st class: 

{xo,Xi} = , {x.i,Xj} = tijkXk, 

which include mass shell condition and requirements of the 
form "spin of the top is equal to the spin of the string" 1 : 

Xo = — Lo = 0, Lo = ^ a* a„, 
X3 = S 3 - A 3 = 0, A 3 = X "' 




0, 



(3) 



X-=S-- A. 
= Si±iS 2 , 



where Si = S e h is a projection of S onto ej. Constraints 
generate the following transformations: 

• Xo generates phase shifts of oscillator variables 

Eq : a n —> o„e~™ T and translations of mean coordinate 
— > Z M + P^r/n; 

• X3 generates phase shifts R 3 : a n —> a n e~ la and rota- 
tions of e*i,2 about e 3 : e% — > R(e 3 , —a)ei, i = 1, 2; 

• Xi,2 generate rotations of basis Si about axes ei,2 and 
certain non-linear transformations of oscillator vari- 
ables, which will be considered in more details in sub- 
sectional 

The main difference of this mechanics from the standard 
light cone gauge description consists in the fact that the 
gauge axis is included in the set of dynamical variables. 
Its position in the center-of-mass frame is described by unit 
norm 3- vector e 3 , while the gauge axis in Minkowski space is 
defined by light-like 4- vector n~ = (N° — N^e l 3 )/V2, where 
N£ is an orthonormal tetrad introduced in Part I. In famil- 
iar terms, Xo generates the evolution of the string, i.e. can 
be used as Hamiltonian; \ 3 transforms the canonical vari- 
ables of the theory in such a way that observable string's 
coordinates and momenta are not changed, therefore it can 



(1) 



1 In this paper the terms "generator of rotations" , "orbital mo- 
ment" and "spin" refer to the same variable. Finally it defines the 
spin of elementary particles, which will result from quantization. 



be interpreted as internal gauge transformation of the the- 
ory; xi,2 change the orientation of the gauge axis S3 and 
are followed by reparametrizations of the string. Generators 
of Lorentz group defined in Part I produce "rigid" Lorentz 
transformations of the world sheet together with the gauge 
axis, without change of its parametrization. As a result, in 
quantum mechanics Lorentz group was freed of anomalies. 
The central problem of this approach is that the algebra of 
constraints \i in quantization acquires the same anomaly 
that earlier was in Lorentz group. To solve this problem, we 
partially eliminate the gauge degrees of freedom in the clas- 
sical theory, imposing additional gauge fixing conditions to 
Xi- These conditions are equivalent to fixation of gauge axis 
e.3 with respect to other dynamical vectors in the system. As 
a result of this procedure, the set of constraints in the theory 
will be reduced, and the remainder will be anomaly-free. 

Before to proceed, we will consider one remarkable par- 
ticular solution in string theory. 

Straight-line string is a solution, for which the string has 
a form of straight segment rotating about its middle at con- 
stant angular velocity in the center-of-mass frame. This so- 
lution plays a special role in string theory 

• it can be quantized at arbitrary number of dimensions 
d > 3 without anomalies in Lorentz group; 

• it belongs to a border of classical Regge-plot P 2 /2iv > 
S, and in quantum theory corresponds to a leading 
Regge-trajectory P 2 /2tv — S; 

• it describes satisfactory the spectrum of light mesons 
and gives an opportunity to consider their radiative 
transitions. 

In light cone gauge the straight-line string is represented 
by one-modal solution, if the gauge axis e$ is directed along 
the spin S, orthogonally to the plane, where the straight-line 
string is rotating: 

On=0, n/1, S± = 0, g = S = S 3 = |ai| 2 . (4) 

There is a gauge-equivalent one-modal solution with exited 
(a_i)-mode and £3 opposite to S. Other directions of gauge 
axis for the straight-line string give infinitely-modal solu- 
tions (particular example see in Part I). Later we will use 
the solution J3J to study the structure of general theory in 
its vicinity. 

Additional gauge fixing conditions we impose have a form 



a a + d-s = 0, a* + a1 s = 



(5) 



for some s > 0. It's easy to check that straight-line string 
solution HI satisfies these conditions at s > 1. Later, in sub- 
section [521 we will show that conditions © can be imposed 
on any solution of string theory. 

The gauge fixing conditions @ are preserved by transfor- 
mations R3. They are not preserved by Eo, however, there 
is a remainder of _Eo- s ymmetry, discrete transformation 



D 2 



t/s 



preserving ©: a s + a_ s — > — a s — a-s- As a result, R3,D 2s - 
symmetries are present in the theory after gauge fixation. 

Substituting © to the symplectic form we see that 
a± 3 -terms cancel each other: -^{da* A da s — da! s A da- s ) = 



0. To study this property in more detail, we introduce the 
variables qi — Re(a s + a_ s )/2, q 2 = — Im(a s + a_ s )/2, pi — 
Im(a s — a_ s )/2, p 2 = Re(a s — a_ s )/2, and rewrite a± s -terms 
of the symplectic form as | (dpi A dqi — dp 2 A dq 2 ) . Gauge 
fixing conditions (J^J are rewritten to q\ — q 2 — 0, now we 
see that the variables pi, 2 , canonically conjugated to qi. 2 , 
drop out from the symplectic form. pi t2 should be expressed 
in terms of other dynamical variables from ^-constraints, 
and independently on how complex these expressions could 
be, the symplectic structure of the mechanics remains to be 
simple. This exclusive property follows from the fact that 
two gauge fixing conditions are in involution with each 
other: {a 3 + a_ 3 , a* + al s } = 0, i.e. generate abelian group 
of transformations. 

Then we see that contribution of a± s to A3 vanishes sim- 
ilarly: i (|a s | 2 — |a- s | 2 ) = 0, as a result, a± s -terms also drop 
out from X3 = S3 — A3. This result follows from the fact that 
gauge fixing conditions @ are preserved by transformation 
R3 and therefore are in involution with X3- 

The gauge fixing conditions @ are not in involution with 

Xo, and a± s -terms in Lo do not vanish: |a s | 2 + | a, s | 2 = 

2|a s | 2 . Poisson brackets of JSJ with x± a l so do not vanish. 
We conclude that © are gauges for (x±> Xo) an d a± s should 
be determined from these three constraints. 



2.2 Algebraic properties of lia-lcg 

Isolating contribution of a± s -oscillators in and using the 
relation a_„ = — a», we have 



a 2 d + |a s a*d* + a s f + a* s g + E_ 




0, (6) 



a* d* + \a s a*d + a*/* + a s g* + E + — \ S-S+ = 



— ' + 2a s a s , 

where 

d = d+-d-, / = /+-/_, g = g~ - g+, 



(7) 



(S) -E' 



a k ak 



d+ = a 2s /s, d- = aL 2s /s, L [ S> = 
3- = ^ j: a,ka s -k, g+ ~ ^ ^afca_ s _fc, 

£- = y ' \o lk a n a" k+n , E+ = S ' j^a k ana k +„, El = E+, 

Here in the sums terms with a *^ and a±l are excluded. 
Introducing denotations 

A = V^> n s = (A 2 - L«)/2, 
k = E_ + \n a d* - XS-, 

we can treat JSJl as overdetermined polynomial system for 
(a 3 ,al) at given values of coefficients (d, f,g,k,n 3 ): 



a 2 d + a s f + a*g + k = 0, 
afd* + a*/* + a s g* + k* = 0, 
a s a* — n s = 0. 



(8) 



To solve this system we find its Groebner basis 0|1]: 

Gr[^,{a s ,a* s }} = {H, a*G + F,a*G + F,...}, (9) 



where H,G,F,... are polynomials of coefficients, given in 
Appendix 1. Groebner basis has the same set of roots as 
original system (SJ and possesses a structure convenient for 
their sequential determination. The first equation H = is a 
condition of consistency of the original system formulated in 
terms of its coefficients. The next equations give the solution 
of the system as rational expressions 

a* =-F/G = -F/G=... 




Fig.l. Resolution of singularities in a*. 

Different expressions for a* s coincide on the surface H — 0, 
in the regions where the denominators do not vanish. In the 
points where the denominator in one formula vanishes, e.g. 
G = 0, we can use the other formula a* = —F/G, provided 
that G / 0. Therefore, in such points both numerator and 
denominator of a* = ~F/G must vanish simultaneously and 
resolution of this ambiguity gives finite value (—F/G). At 
the end of Groebner basis the polynomials of degree > I 
in a* are present. Those which do not contain a s , describe 
solutions for the exceptional degenerate cases when all de- 
nominators of the linear formulae vanish: G = G — ... = 0. 
Those which are linear in a s formally can be used to find 
a s . The original equations are also appended at the end of 
Groebner basis, meaning that all solutions must be checked 
by substitution to the original system. 

The presence of complex conjugation in JSJ creates two 
general problems: (i) the described procedure treats (a s ,at) 
as independent variables (a B ,b B ) £ C 2 , it does not guar- 
antee that a s will be complex conjugated to b s . However, 
the symmetry of original system with respect to the com- 
plex conjugation implies that (b*,a*) is also a solution of 
@, possibly the other one. We are interested only by the 
case (6s, a^) = (a s ,b s ). (ii) theorems |3] in the case of non- 
holomorphic polynomial systems do not guarantee that all 
solutions of H = will be continued to full solutions of JSj, 
obstacles are the equations of type z* z — — 1. To resolve 
both problems we substitute the expression al = —F/G and 
its complex conjugate a 3 = —F*/G* to the original equa- 
tions and explicitly check that they are satisfied on the sur- 
face H = (after putting all terms over a common denomi- 
nator we check that H is a polynomial divisor of the numer- 
ator). As a result, we obtain the solutions (a s ,a* s ), which 
are correctly transformed by complex conjugation and are 
defined everywhere except of zero-measure set of degenerate 
cases described above. 

H = 0, the condition of consistency of the system 10, is 
the constraint, remaining after imposition of two gauge fix- 
ing conditions 10 to three initial constraints xo,X±- This 
condition is real-valued: H* — H. Being used as Hamilto- 
nian, it generates correct string evolution, consisting of shifts 
a — > a + t and such reparametrizations that keep gauge fix- 
ing conditions permanently satisfied. Hamiltonian is in 
involution with the constraint X3, i.e. is preserved by trans- 
formation 7?3 : a n — > a n e~ ta , S- —> S-e~ la . Hamiltonian is 



also preserved by transformation Dis : a„ — > o„e~™' s . Us- 
ing the explicit representation of a s from Appendix 1, one 
can easily verify that the properties of a s with respect to 
these transformations are the same as before gauge fixation: 

i?3(a s ) = a s e~' a , _D 2s (a s ) = -a s . (10) 

It's important to remark that these symmetries are defined 
by polynomial structure of the expressions: the symmetries 
implemented as phase rotations of oscillator variables im- 
pose certain conditions on the powers of oscillator variables 
in monomials. As a result, R3, 7?2s-symmetries can be pre- 
served in quantum theory as well. 

Generally H is 8th degree polynomial in A, in special 
case S± = it degenerates to 4th degree polynomial in A 2 . 
Using implicit function theorems, one can locally rewrite 
the constraint H = in the form of mass-shell condition: 
P 2 — $ = 0, where the function $ globally has multiple 
branches. Because at s > 1 straight-line string J1J satisfies 
gauge fixing conditions one branch of $ passes through 
the straight-line string. Further we study in details the struc- 
ture of solutions in the vicinity of straight-line string. 

2.3 Solutions in the vicinity of straight-line string 

Lemma 1: at s — 2 in the vicinity of straight- line string the 
system © has unique solution, representable as C°°-smooth 
analytical function of coefficients of the system. 
[Proof of the lemma can be found in Appendix 2.] 

Let's fix s = 2. On the straight-line string solution |gj for 
the coefficients of the system © we have f± = d± = S± = 
S± — 0, in g the term g+ = 0, while g~ = a 2 . Let's pull off 
ai-contribution from g-terms: gL = g~ — a 2 , g — g'_ — g+. 
Due to the lemma above we can represent 02 as Taylor's 
expansion with respect to these coefficients. Further we will 
study a part of this expansion, which do not contain Si- 
terms. The case S± = possesses convenient simplifications, 
particularly, the first two equations in (|HJ define a closed 
system for (02,02) and decouple from the third one, which 
gives direct definition of P 2 : 

£ = 4 2) +2|a 2 | 2 . 

Lemma 2: (S± = 0)-terms of 02 expansion in the vicinity 
of straight-line string have a form 

n>l 

where P n are polynomials of 01, S_, d, f, g' and their con- 
jugates. The polynomials possess the following properties: 
transformation D4 : a n — > a n i~ n preserves the monomials of 
P n , transformation Rg : a n —* a n e~ la multiplies the mono- 
mials to e la ; each monomial in P n contains, counting powers, 
n variables from the group (E_ , d, f, g' and their conjugates) 
and 2(n — 1) variables from the group (ai,^); each mono- 
mial contains at least one variable E± ; numbers of variables 
from (£±)- and (d, d*)-groups are related: n(E) = 71(d) + 1. 

Explicit expressions for the first six polynomials P n can 
be found in Appendix 1. The expansion il l U has an obvious 
singularity at ai = 0. Later we will show that in quantum 
theory this singularity can be removed by a particular choice 
of quantum ordering. 



2.4 Geometric properties of lia-lcg 

In this subsection we describe in more details the transforma- 
tions, generated by constraints Xo*Xi- At first, we introduce 
several definitions. 

Let's consider in 3D space: smooth closed curve Q(u) with 
marked point O and unit norm vector S3. Let's introduce 
variables 

. exp [2jj» (L(a) - (Q(a) - Q(0))e 3 )] , (12) 

where I/(<r) is a length of arc between points O and Q(o) 
along the curve, L t is total length of the curve. Two proper- 
ties obviously follow from the definition: a_ n = So — 0. 
Let's decompose vectors a n into the components, parallel 
and orthogonal to S3: a n = an3e3 + a n x, and denote real and 
imaginary parts of S n± as q n± and p nA _ : S n± = q n± + ip n ± ■ 
Let's fix some n — s > and write g s x = q and p s x = P- 
Functions q(S-j,),p(S3) define smooth vector fields on unit 
sphere of S3 (tangent to the sphere). Due to topological 
"hedgehog" theorem, these fields have singular points on the 
sphere, where corresponding field vanishes, e.g. q — 0. 




Fig. 2. Definition of vector fields q(e3), p(e3). 

Note: the curve Q is a projection of supporting curve, in- 
troduced in Part I, into 3-dimensional subspace, orthogonal 
to total momentum P M . In terms of string dynamics, it is a 
trajectory of string's end in the center-of-mass frame. The 
length of this curve is equal to double mass of the string: 
Lt = 2V P 2 . Gauge axis S3 relates the following parametriza- 
tion to this curve: 

2ir 

a= — (L(a)-Q 3 (a) + Q 3 (0)), (13) 

where Q3 = QS3. Now we recognize in 1121 Fourier modes 
of function a(a) — Q'(a), where a is leg-parameter 1131 . 
This expression is written in parametric-invariant form, as 
circulation integral. Vectors a n ± are related with earlier 
introduced oscillator variables follows: 

a n ± — a n \S\ + a n iSi, 

a n i = (o n + al n )/2, a n2 = i(a n ~ a*„)/2. 

Lemma 3: X3 generates transformations, preserving S3 and 
OnJ_; Xi,2 generate rotations of S3 and associated changes 
of a n ± according to the formula 1121 . Gauge fixing con- 
ditions (Sj correspond to a singular point q = 0. Evo- 
lution, generated by xo, is represented as phase rotations 
a n ± — > a n ±e~ lnT . During this evolution vector uj — q x p is 
preserved, and singular points q = move along zero-level 
curves of a function F(e*3) = LJS3 = 0. 

This lemma implies that gauge fixing conditions (JKJ can 
be imposed on any solution of string theory. Namely, these 



conditions can be satisfied for any curve Q(o~), directing the 
light cone gauge axis S3 to a singular point qXSs) — of a 
vector field on the sphere, constructed in terms of this curve. 

In general position the vector fields on the sphere have 
even number of singular points, which is > 2. Non- 
degenerate singular points of 2D vector fields are 




center focus node saddle 

Fig. 3. Singular points of 2D vector field. 



For each type index of singularity is defined as algebraic 
number of rotations of vector <f(e*3), when S3 passes around 
singular point (> 0, if directions of rotations of q and S3 co- 
incide; < otherwise). For (center, node, focus) it is +1, for 
saddle —1. The sum of all indices is equal to Euler charac- 
teristic of the surface V — E + F: (num. of vertices) — (num. 
of edges) + (num. of faces) for any tessellation of the sur- 
face into faces, which for the sphere equals 2. Thus, generic 
vector field on a sphere has 2 singular points of index +1 
and arbitrary number of self-compensating pairs (+1,-1). 

Presence of several singular points q\Sz) = leads to 
the fact that an orbit of a gauge group generated by \- 
constraints intersects a surface of gauge fixing conditions 
© in several points of the phase space. This phenomenon 
is also encountered in the theory of non-abelian gauge fields, 
where it has been studied by V.N.Gribov 0- We will call 
such points Gribov's copies. 




Fig. 4. Gribov's copies: the orbit of a gauge group o(G) 
intersects the surface of gauge condition 5 = in several points. 

In our case the copies comprise different sets of oscilla- 
tor variables {a n }, reproducing the same curve Q(cr) (differ- 
ently parameterized), i.e. physically correspond to the same 
string's motion. This property indicates remaining discrete 
gauge equivalence in the theory. To identify the equivalent 
states, in classical theory the phase space should be factor- 
ized with respect to this symmetry. Analogous procedures 
can be applied in quantum mechanics, e.g. by constructing 
irreducible representations for operators possessing this dis- 
crete symmetry and formulating respective selection rules. 

During the evolution Gribov's copies move along closed 
loops, constructed by lifting of zero-level curves F{Sz) = 
to the space of oscillator variables. In more details, this 
construction is performed as follows: vector e*3(r) moves 
along the contour F{ez) = on the sphere, defined for 
a given value s, while for n 7^ s expression (II 2t defines 
time-dependent vector fields a n x(e3, t), where we substitute 
S3 — S3(t). Note that corresponding vector fields 
generally do not vanish on these loops. The set of Gribov's 
copies possesses D2s-symmetry, however in general case this 
symmetry does not exhaust the discrete gauge equivalence. 

Lemma 4: during the period of evolution At = 2-7T singular 
points q = return to the initial position 2s times. The 
lifting of the contour F(Ss) = to the space of oscillator 



variables has at least 2s Gribov's copies, related by D23- 
symmetry. 

Below we study the structure of Gribov's copies in the 
vicinity of straight-line string solution J3J and its S3 — * —S3 
companion. For these purpose we investigate the expan- 
sion of q(e3) for this solution in the vicinity of northern and 
southern poles of the sphere S3 = ±S/\S\. 

Lemma 5: <j(e3) for straight-line string solution at s = 1 
has no singularities in the vicinity of the poles; at s > 1 
has in each pole a singularity, which does not move during 
the evolution, and at s — 2 is a saddle point, at s > 2 has 
degenerate type. 

Lemma 6: Solutions, close to Q, have at s = 2 a sad- 
dle point moving along a small loop in the vicinity of the 
northern pole. During the period of evolution At = 2n the 
saddle point performs 4 turns. Associated lifting has exactly 
4 Gribov's copies, related by D4-symmetry. 

Analogous structure of Gribov's copies is established for 
the companion solution near the southern pole. Solutions 
close to straight-line string at s — 2 should also have other 
singularities faraway of the poles (saddle points have index 
— 2, while the total index of singular points on the sphere 
should be +2). Further in quantum mechanics we will use 
the expansion 1111 in the vicinity of solution JIJ, correspond- 
ing to the northern pole singularity. 

Summary of classical mechanics: the phase space 
(Z, P, a n , a*n i S, Si) possesses trivial topology. Hamiltonian 
H is globally defined in it as a polynomial function (Ap- 
pendix 1). The following two features can be considered 
as topological defects: (i) the expression of mass squared 
P 2 — $ from the constraint H — is multi-valued; (ii) there 
is a discrete remainder of original gauge invariance, which 
includes 7?2s-symmetry, but generally is not exhausted by 
it. In the vicinity of straight- line string solution @ at s = 2 
the branch of P 2 = $ is regular, and the discrete gauge 
invariance is exhausted by D4-symmetry. 



3 Quantum mechanics 
3.1 Canonical operators 



iw>= n 1 
n 



K4) 



fc<0,fc^-s 



k>0,k^s 
1 



K)" fc |o)- 



For instance, we will write |li2_i^ = ^=afa?_i|0^ etc. 
So defined state vectors have positive norm. Occupation 
numbers 



rik = 



k\ - akak -- \k\ 



a^at, k > 
auat, k < 



= 0,1,2... 



Note: one-component oscillator variables a k used through- 
out this paper are related with commonly applied two- 
component oscillators dki,2 by the formulae of subsec- 
tion 1231 which on quantum level become: at = dki — ia,k2, 
ajj" = a_fci + ia_fc2- The inverse formulae are aui = 
{a,k + oi fe )/2, at2 = iiflk — a ^fc)/2. Here taking two sets of 
oscillator variables Ofc», i = 1,2 with at = a-ki we construct 
one set without this property. Usage of one-component 
oscillators simplifies the algebra. It's easy to verify that 
definition of vacuum 1151 is equivalent to a standard one 
«fci|0^) = 0, k > and the states are the linear com- 

binations of ILx^m)"" (4 2 )" fc2 10>. 

Quantum top: the space of states is formed by functions 
^(e), e\ G 5*0(3). Because the rotation group is not simply 
connected, 7Ti(SO(3)) = Z2, two representations are possi- 
ble: single- and double- valued !7.. Spin is defined as dif- 
ferential operator S l = — ieyfeef d/de* , while the projection 
of spin onto the coordinate system e*; is Si = ieijkeid / 'de\. 
Operator S l generates the rotation of the coordinate system 
Si in external space, while Si generates the rotations about 
the axes ei . These transformations act on different indices in 
ej. They commute, therefore S 3 and 53 are simultaneously 
observable. 

The basis in the space of states is formed by Wigner's 
©-functions 



|55 3 5 3 ) =2>i 3S3 (e), S 3 ,S 3 



-5,-5- 



[Zfj,,Pjy] — —igpv, 

[a k , at] = k5kn, fc,n#0,±s, (14) 
[S\S J ] = ie ijk S k , [S\e)] = ie ik ie l ], 
[Si,Sj] = —itijkSk, [Si,e k ] = —itijiei, 
[5\5,]=0, efej = *j, S^elS 3 . 

The space of states is a direct product of three components: 



where 5 are integer for single- valued representation and half- 
integer for double-valued one. Further a constraint S3 = 
A3 £ Z will select only integer 5- values 2 . 

Operator 5 2 = 5*5 l = SiS t has eigenvalues 5(5 + 1). 
The constraints, being functions of Si and oscillator vari- 
ables, commute with S 2 , as a result, the determination of 
mass spectrum can be performed separately for each 5 value. 
Particularly, spin singlet 5 = corresponds to state vectors 
independent on e\. On such states all 5 l = Si = 0. 



Space of functions \P(P) with the definition Z — —id/dP. 
Fock space with a vacuum 

o fc |0} = 0, fc > 0, a£|o) = 0, fc<0 (15) 
and states created from vacuum by operators 



Constraints X3|^) = ^ p^) ~ are imposed on the states 
in the spirit of Dirac's quantization of constrained theories. 
Here X 3 = 5 3 -4 S) > 

2 However, using formal modifications A3 — > A3 + 1/2, classi- 
cally vanishing after the recovery of Planck's constant, it's possi- 
ble to obtain half-integer spin as well. The question on possibility 
of inclusion of half-integer S-values in this mechanics should be 
better investigated. 



A 



^2 % ■ a t a * ■■ = E sign k 



fc^O,±s 



In the next subsection we will constr uct two versions of 
Hamiltonian H: polynomial of A = * / P 2 /2ir and P 2 — $ 
with $ is defined by Taylor's series li lt . For the opera- 
tors entering in H: d±, f±, p±,£± and their conjugates we 
reserve a term elementary operators. They can be defined 
by their classical expressions with only replacement of 
complex conjugation to Hermitian one. These definitions 
have no ordering ambiguities. For the states with the finite 
number of filled modes the matrix elements of elementary 
operators are described by finite sums. 



Symmetries R3, Eo, Dis in quantum theory are generated 
by operators 

R3 = e Xi , Eo = e , D 2s = e ' . 

It's convenient to introduce the notion of (AA3,ALq)- 
charges for operators satisfying the relations: 



[4 S) , op] = AA S ■ op, [L 31 , op] = ALq ■ op 



i.e. the operators, which increase or decrease the quantum 



numbers (A ( 3 a \L 3 ') by (AA 3 , AL ) units: 



00 



op 



4 s) ,4 s) ) = |4 s) +AA 3 ,4 s) +ALo). 



Particularly, AA3(a n ) — —1, ALo(a n ) = —n. Hermitian 
conjugation of operators reverses their charges. The el- 
ementary operators have definite (AA3, ALo)-charges, see 
Table 1. The discrete symmetry D2 S is also characterized by 
the charge Q = Lq S ' mod 2s, so that the subspaces with given 
Q are eigenspaces Z?2s|Q) = e l7r ®/ 3 \Q}, and D2s-symmetric 
operators keep these subspaces invariant. 

Because the symmetries of classical theory in our case 
are defined by its polynomial structure, they are preserved 
on the quantum level. Particularly, quantum Hamilto- 
nian possesses -R3, Z^s-symmetries, as a result, its non-zero 
matrix elements form blocks, located on Ag^Q-diagonal: 
(A^Q\H\A ( ^Q). 



3.2 Quantum Hamiltonian 

Hamiltonian H can be defined directly by the expressions 
of Appendix 1, where we only need to fix a proper ordering 
(Table I). As a result, we obtain a polynomial spectral prob- 
lem of the form H(X)^ X = 0, where H(X) = £ h n X n , the 

coefficients h n are Hermitian operators and A = y/ P 2 /2ir is 
a spectral parameter. The problems of these kind are equiv- 
alent to a standard linear spectral problem in larger space: 
(M - A)* A = 0, where 
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I is a degree of polynomial H(X) and h n = h l 1 h n , pro- 
viding that exists. Using the regularizations based 
on the restriction of operators to finite-dimensional spaces 
|Lq S ' < (described below), the problem is reduced to a 
spectral analysis of a large matrix, which can be performed 
with the aid of standard computational algorithms Due 
to the block structure of Hamiltonian, A spectrum can be 
defined separately in each A3 Q-block. 

The following circumstance is rarely remarked. Generally 
Hamiltonian H can depend on P 2 non-linearly. Classically 
one can resolve the constraint H — with respect to P 2 , 
and obtain a locally equivalent theory. On the quantum 
level these theories have principial differences: the first one 
creates non-linear spectral problem, while the second one 
gives the spectral problem (P 2 — ^)'i> = 0, linear in terms of 
P 2 . The last problem has relatively simple structure of solu- 
tions: Hermitian operator $ possesses complete orthonormal 
basis of eigenvectors and real eigenvalues, representing the 
mass spectrum. In the polynomial problem Hermitian oper- 
ator H(X) at fixed A also has real eigenvalues and complete 
orthonormal basis of eigenvectors: H(X)'i/\.f 1 = ^(X)^ \ :fl . 
Here the eigenvalues and eigenvectors are the functions of A. 
However, the vectors satisfying H(X)ty\ — 0, defined by the 
roots of the equation /x(A) = 0, should not necessarily form 
a complete basis. For finite-dimensional matrix systems the 
number of real roots of this equation can be less than the 
dimension of the space, it also can be greater than the di- 
mension. Furthermore, the vectors = for different A 
generally are not orthogonal. 

The last property can be improved. In a particular case 
when for the real roots A form a complete basis, we can 
construct the matrix Y — AAA -1 , where the columns of 
X are vectors ^a and A = diag A. It's easy to verify that 
under these conditions Y is a solution of matrix equation 
h n Y n = 0, where Y n stand on the right of h n , while Y + 
is a similar solution from the left: ^2Y +n h n = 0. In the 
classical limit both Y and Y + correspond to a real root A 
of classical equation ^/t n A n = 0. Operator (Y + Y + )/2 is 
Hermitian quantum analog of A. It has orthogonal system of 
eigenvectors and can be used to define the mass spectrum. 

The problem with the number of eigenvalues is more seri- 
ous. If the classical equation H(X) — has several real roots 
Xi, one can represent each root by a separate function on the 
phase space. In quantum mechanics each root becomes Her- 
mitian operator possessing a complete system of eigenvec- 
tors. However, in the solutions of polynomial spectral prob- 
lem the corresponding eigensystems are mixed together and 
their separation appears to be non-trivial problem. Theory 
of matrix polynomials factorization |Sj possesses sufficient 
potential to solve this problem. Particularly, it introduces 
sign characteristics of eigenvalues, which under certain con- 
ditions allow to separate them completely. On the other 
hand, the resolution of constraint H = to P 2 — $ and 
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Table 1: L Q -normally ordered elementary operators 
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separation of $ branches on the classical level provides an 
alternative approach to this problem. Particularly, the spec- 
trum of P 2 in the vicinity of straight-line string can be stud- 
ied using the expansion 1111 , In this paper we will consider 
this problem for a particular case of spin singlet S = 0. 

3.3 Solutions in the vicinity of straight-line string 

Let's fix s — 2 and put S3 = S± = in the formulae above. 
Quantum analog of (II li t is constructed as follows: 



2 



E 



D~ 



(16) 



where D = afai + 1/2 = (afai + aiaf)/2. 

Polynomials P n are defined by expressions of Appendix 1, 
where we substitute the definitions @ of elementary oper- 
ators and fix the ordering, shown in Table 1. This ordering 

puts Lq -lowering elementary operators to the right from 

(2) 

L -raising ones, thus providing better convergence proper- 
ties for the expansion (1161 . Particularly, the matrix elements 

of a,2 between the states with finite L 2 ' are given by finite 
sums, and large n terms of 1161 contribute only to the matrix 

(2) 

elements with large L , due to the following lemma. 
Lemma 7: 

(L^ =N 1 \:P n :\tf ) 



= N 3 ) = 0, if 4(n - 1) > Ni + N 2 . 



Square-mass operator is given by 



— = L (2) +2aja 2 . 

ATT 



(17) 



One can also include here classically vanishing constant 
term: P 2 — » P 2 +c. In this definition (at c > 0) the operator 
P 2 is Hermitian and positively defined 3 . 

The introduced ordering possesses the other feature conve- 
nient for computations: the matrix elements of the normally 
ordered operators restricted to finite-dimensional subspaces, 
behave regularly when the dimension of subspaces increases. 
Namely, the matrix in larger space includes the matrix in 
smaller space as an exact submatrix. Without normal or- 
dering this simple behavior would be violated. 

Lemma 8: let opi(N) = (L { 2) < N\opi\L {2) < N) be re- 
strictions of elementary operators opi to finite-dimensional 



space Lq < N. Let : P(opi) : be L z; -normally or- 
dered polynomial of opi. Then (i) at Ni < N2 the ma- 
trix : P(opi(Ni)) : is a submatrix of : P(opi(7V 2 )) : and (ii) 
: P(o Pl (N)) := (L (2) < N\ : P(o Pl ) : \L {2) < N). 

Practically, we compute the matrix elements of elemen- 

(2) 

tary operators in L < N subspace. The dimension of this 
subspace is rapidly increasing with N, e.g. dim — 106587 



'(2) 



3 Note that other possible definition P 2 /2n ■ 
does not have this property. 



•(2) 



a, a2 



for iV = 25, so that the elementary operators are rep- 
resented as matrices of very large size (dim x dim), see 
Appendix 3. The matrices have noticeable block struc- 
ture, corresponding to their (A^ , L 2 ' )-charge properties. 
As a result, non-zero matrix elements of 02, necessary for 
computation of P 2 at S = 0, are located in the blocks 
(yl 3 = -l,Q + 2|a 2 |-43 = 0,Q), while P 2 itself is located in 
(^Az — 0,Q|P 2 |j43 = 0,Q^. In addition to these properties 

we use the fact that elementary matrices inside (A3 2 ' , Z/ 2 ' )- 
blocks are very sparse (at large their non-zero content 
occupies less than 1% of the blocks), and implement special 
algorithms for sparse block matrix computations, described 
in more details in Appendix 3. Finally, we determine the 
spectrum of P 2 /2vr up to the value N = 28 (dim = 260256) 
and the number of terms in expansion 1161 up to n — 6. The 
initial part of P 2 spectrum is rapidly stabilized when n is 
increasing. Particularly, the region P 2 /2n < 10 is defined 
by n < 3 terms, while n > 3 corrections influence the higher 
P 2 values only. Figure |H| shows the resulting spectrum of 
P 2 , for separate Q-values and superimposed. 
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Fig.5. Spectrum (P 2 /2n, S = 0). 
Table 2: eigenvectors (P 2 /2tt < 5,S = 0). 
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The spectrum contains degenerate eigenvalues, whose 
multiplicities are shown on the figure by numbers in brackets 
(non-degenerate eigenvalues have no these numbers). Black 

points display the states annulated by operator a 2 , which are 

(2) 

simultaneously the eigenvectors of L and therefore have 
integer P 2 /2tv. The beginning of the spectrum consists only 
of such states, as the following lemma shows. 



Lemma 9: a 2 annulates \L, 



(2) 



<5). 



At P 2 /2-7T > 6 the eigenvalues of L 2 ^ start to split. Par- 
ticularly, one state with L Q 2> = 6 goes to P 2 /2tt = 8.16792, 
two states with L {2) = 8 become P 2 /2vr = 9.0107, 10.018 
etc. Such states (displayed by white points on figEJ cor- 
respond to the exited a 2 -modes. The values of P /2ty for 



these states are not integer (however, some of them are nu- 
merically close to integer). Many eigenvectors of > 6 
are also annulated by 112, their P 2 /2ir are degenerate and 
integer. This degeneracy can be removed by further model 
corrections, such as spin-orbital interaction discussed in 

Subspaces with different Q in quantum theory play a role 
of factor-spaces with respect to D4-symmetry, which is sat- 
isfied exactly for Q — 0: = and in projective sense 
for Q / 0: ~ ty. Quantum Hamiltonians of interac- 
tion, possessing D4-symmetry, keep these subspaces invari- 
ant. Hamiltonians, violating this symmetry, create transi- 
tions between different Q-sectors. If the violation of D4- 
symmetry is small, the particles on the first daughter Regge 
trajectories will be quasi-stable. 



[11] W.H. Press, S.A. Teukolsky, W.T. Vetterling, B.P. 
Flannery, Numerical Recipes in C, Cambridge Univer- 
sity Press, 1988. 
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Summary of quantum mechanics: representation of canon- 
ical operators 1141 is straightforward. Quantum Hamil- 
tonian can be defined either as polynomial function (Ap- 
pendix 1) or as expansion Ijlfip in the vicinity of straight-line 
string, both properly ordered. In the first definition a non- 
linear spectral problem appears, with associated separation 
task for the roots of characteristic equation. The second def- 
inition leads to the ordinary linear spectral problem, which 
allows, particularly, to determine the mass spectrum for the 
states of spin singlet 5 = 0. 

A paper devoted to systematic study of Gribov's copies in 
this mechanics and determination of mass spectrum outside 
of the spin singlet will be published. 
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Appendix 1: the polynomials, entering to the definition of Groebner's basis © of the system JSJ: 



H = -/* gg* 2 k 2 -fg 2 g* k* 2 - d f* g k k* 2 - d* f g* k 2 k* + d* g* 2 k z + d g 2 k* 3 + d* f f*g*kn s +dff* 2 gk*n a + 2d* f 2 gg* k* n 8 + 
2df* 2 gg*kn s -3d*fgg* 2 kn s -3df*g 2 g*k*n s -2dd*f*g*k 2 n s -2dd*fgk* 2 n s -d 2 f* 3 gn 2 -d* 2 f 3 g* n 2 s +d 2 d* f* 2 kn 2 s + 
dd* 2 f 2 k* n 2 + 3dd* 2 f g* kn 2 + 3d 2 d* f* gk* n 2 + f f* g g* kk* +g 2 g* 2 kk* + dd* k 2 k* 2 - g 3 g* 3 n a +2 f f* g 2 g* 2 n s - f 2 f* 2 g g* n s - 
dd* f f* k k* n s + dd* g g* k k* n a - dd* f f* g g* n 2 + 3 dd* g 2 g* 2 n 2 - 2d 2 d* 2 kk* n 2 - 3d 2 d* 2 g g* n 3 - d 2 d* 2 f f* n.| + d 3 d* 3 n 4 

F = -f* 3 g 2 g* kn s + 2d* f* 2 g g* k 2 n s - d* 2 f* g* k 3 n s + f* 2 g 3 g* k* n s - d* 2 gg*k 2 k*n s - d* g 3 g* k* 2 n s + df* 4 g 2 n 2 a - 
d* f*g 3 g* 2 n 2 - 2dd* f* 3 gkn 2 - 2d* 2 f f* g g* kn 2 + d* 2 g 2 g* 2 kn 2 + dd* 2 f* 2 k 2 n 2 + 2d* 3 f g* k 2 n 2 - 3dd* f* 2 g 2 k* n 2 + 
2d* 2 f g 2 g* k* n 2 + Add* 2 f* g k k* n 2 - dd* 3 k 2 k*n 2 + dd* 2 g 2 k* 2 n 2 + 2dd* 2 f f* 2 gn 3 - d* 3 f 2 g g* n 3 s + 2 dd* 2 f* g 2 g* n 3 - 
2dd* 3 f f* kn 3 -2dd* 3 g g* kn 3 a -2dd* 3 f gk* nf + dd* 4 f 2 n 4 - d 2 d* 3 f* gn 4 + d 2 d* 4 knj 

G = f* 3 g 2 kk* - 2d* f* 2 gk 2 k* + d* 2 f* k 3 k* - f* 2 g 3 k* 2 + d* 2 gk 2 k* 2 + d* g 3 k* 3 - f f* 4 g 2 n s + f* 3 g 3 g* n s + 2d* f f* 3 g kn s - 
3d* f* 2 g 2 g* kn s - d* 2 f f* 2 k 2 n s + 3 d* 2 f* g g* k 2 n s - d* 3 g*k 3 n s + 3 d* f f* 2 g 2 k* n a - 2d* 2 f f* gkk* n a - d* 3 f k 2 k* n s - 
3d* 2 fg 2 k* 2 n s - 2 d* 2 f 2 f* 2 g n 2 + 2 d* 3 f 2 f* kn 2 + 3 d* 3 f 2 g k* n 2 - d* 4 f 3 n 3 

The polynomials, entering to the expansion (I1H of variable ci2|s ± =o in the vicinity of straight-line string solution P): 
Pi = -E+, 

P 2 = a* 2 /*E_ +ai 2 g'*S+, 

2P 3 = -2a* 4 /*g'E_ -2a 1 2 a* 2 f* g'* E_ - 2 af d* E_ 2 - 2 ai 2 a* 2 f f* E + - 2ai 4 g'* 2 E+ - ai 2 a* 2 dE_ £+, 

2P 4 = 2ai 2 a* 4 //* 2 E_ + 2af f*g' 2 T,- + 2 ai 2 a* 4 /* g 1 g'* E_ + 2 ai 4 a* 2 f* g'* 2 E_ + ai 2 a* 4 df* E_ 2 + Aa* 6 d* g'T,- 2 + 
2ai 2 a* 4 d* g'* E_ 2 + 2 a x 2 a* 4 / /* g> E+ + 4 ai 4 a* 2 f f* g'* E + + 2a 1 6 g'* 3 Y, + + 4 a ± 2 a* 4 d* f £_ S + + ax 2 a* 4 d* /* E_ S + + 
ai 2 a* 4 d 9 'S_ S+ + 2 ai 4 a* 2 d g'* S_ S+ +ai 4 a* 2 d/S+ 2 +2ai 4 a* 2 d/* S+ 2 , 

4P 5 = -Sa! 2 ^ 6 //* 2 ^^. - 4a* 8 /V 3 £_ " 8 a ± 4 a* 4 / /* 2 <?'* S_ - 4a x 2 af /* 9 ' 2 g 1 * S_ - 4 a x 4 a* 4 /* g' 9 '* 2 S_ - 
4ai 6 a* 2 /* g'* 3 S_-12ai 2 a* 6 d*//*2- 2 -2ai 2 a* 6 d* /* 2 S_ 2 - 4 ai 2 a* 6 d /* g' S_ 2 - 12 a* 8 d* g' 2 S_ 2 - 4ai 4 a* 4 d f* g'* S_ 2 - 
8ai 2 a* 6 d* g' g'* S_ 2 -4 ai 4 a* 4 d* g'* 2 S_ 2 -2 ai 2 a* 6 d d* E_ 3 -4 ai 4 a* 4 f 2 f* 2 S+ -4 ai 2 a* a f f* g' 2 S+ -8 ai 4 a* 4 / /* g 1 g'*S+- 
12ai 6 aJ; 2 //*g'* 2 S + - 4ai 8 g'* 4 S + - 8 ai 4 aj 4 d / /* £_ S+ - 8 ai 4 a* 4 d /* 2 S_ S + - 16 a x 2 a* 6 d* / g' S_ S + - 
4ai 2 a* 6 d* /* g'S_ £+ - 2ai 2 a* 6 dg' 2 E_ S+ - 16 ai 4 a* 4 d* / g 1 * E_ S+ - 4ai 4 a* 4 d* /* g'* S_ E+ - 4ai 4 a* 4 dg' g'* £_ S+ - 
6ai 6 a* 2 dg'* 2 S_ £+ -ai 4 a* 4 d 2 S_ 2 S+ - 4ai 2 a* 6 d* 2 S_ 2 S+ -4ai 4 a* 4 d* f 2 S+ 2 - 2 ai 4 a* 4 d* / /* S+ 2 - 2 ai 4 a* 4 d / g' S+ 2 - 
4ai 4 a* 4 d/* g'S+ 2 - 6ai 6 a* 2 d/g'* S+ 2 - 12ai 6 a* 2 d/* g'* S + 2 - 9ai 4 a* 4 dd* S_ S+ 2 - 2ai 6 a* 2 d 2 S+ 3 , 

4P 6 = 4a! 4 a* 6 / 2 /* 3 S_ + 12 a x 2 a* 8 / /* 2 g' 2 E_ + 4a* 10 /* g' 4 £_ + 16 ai 4 a* 6 / f* 2 g> g>* E_ + 4 ai 2 a* 8 /* g' 3 g'* S_ + 
12ai 6 a* 4 //* 2 g'* 2 S_ + 4ai 4 a* 6 /* g' 2 g'* 2 S_ + 4ai 6 a* 4 /* g' g'* 3 S_ + 4ai 8 a* 2 /* g'* 4 S_ + 6 ai 4 a* 6 d / f* 2 E_ 2 + 
4ai 4 a* 6 d/* 3 E_ 2 + 36 ai 2 a* 8 d* / /* g> E_ 2 + 6 a x 2 a* 8 d* /* 2 g' E_ 2 + 6 a x 2 a* 8 d /* g' 2 E_ 2 + 16 a* 10 d* g' 3 E_ 2 + 
24a 1 4 a* 6 d*//*g'*E_ 2 + 4 ai 4 a* 6 d* /* 2 g'* E_ 2 + 8 ai 4 a* 6 d/* g' g" E_ 2 + 12 ai 2 a* 8 d* g' 2 g'* E_ 2 + 6 ai 6 a* 4 d /* g'* 2 E_ 2 + 
8ai 4 a* 6 d* g' g'* 2 E_ 2 +4ai 6 a* 4 d* g'* 3 E_ 2 + 8 ai 2 a* 8 d* 2 / E_ 3 + ai 4 a* 6 d 2 /* E_ 3 + 6 ai 2 a* 8 d* 2 /* E_ 3 + 6 ai 2 a* 8 d d* g'E_ 3 + 
4a! 4 a* 6 dd*g'*E_ 3 + 8 ai 4 af f 2 f* 2 g> E+ + 4 Ql 2 af f f* g' 3 S + + 12 aj 6 a* 4 / 2 f* 2 g>* E+ + 8 a x 4 a* 6 / /* g' 2 g" S+ + 
12a 1 6 a* 4 //*g'g'* 2 E+ + 16 af f f* g'* 3 S+ + 4a 1 10 g'* 5 E+ + 24 aj 4 af d* f 2 f* E_ E+ + 8 a x 4 af d* f f* 2 E_ E + + 
16ai 4 a* 6 df f* g' E_ E + + 16 ai 4 af df* 2 g' E_ E+ + 24 ai 2 a* 8 d* / g' 2 E_ E+ + 6 ai 2 a* 8 d* /* g' 2 E_ E + +2ai 2 a* 8 dg' 3 E_ E+ + 
24ai 6 a* 4 d//*g'*E_ E+ + 24 ai 6 af d f* 2 g" E_ E+ + 32 ai 4 a* 6 d* f g' g'* E_ E+ + 8 ai 4 a* 6 d* f* g' g" E_ E+ + 
4ai 4 af dg' 2 g'* E_ E+ + 24ai 6 af d* f g'* 2 E_ E++6ai 6 a* 4 d* /*g'* 2 E_ E++6ai 6 a* 4 dg' g'* 2 E_ E++8ai 8 af dg'* 3 E_ E+ + 
12ai 4 a* 6 dd*/S- 2 E+ + 28ai 4 a* 6 dd* /*E_ 2 E+ + 2 ai 4 a* 6 d 2 g' E_ 2 E+ + 12 ai 2 a* 8 d* 2 g' E_ 2 E+ + 3 ai 6 af d 2 g'* E_ 2 E+ + 
8a! 4 a* 6 d* 2 g'* E_ 2 E+ + 6 ai 6 af d f 2 f* E+ 2 + 12 ai 6 a* 4 d / f* 2 E+ 2 + 8 ai 4 a* 6 d* / 2 g' E + 2 + 4 ai 4 a* 6 d* f f* g' E + 2 + 
2ai 4 a* 6 d/g' 2 E+ 2 + 4 a x 4 af d f* g' 2 E+ 2 + 12 ai 6 a* 4 d* f 2 g'* E + 2 + 6 ai 6 a* 4 d* / /* g" E + 2 + 6 ai 6 af d f g> g" E+ 2 + 
12ai 6 a* 4 d/*g'g'*E + 2 + 12 ai 8 af d f g'* 2 E+ 2 + 24 aj 8 a* 2 d /* g'* 2 E+ 2 + 3 ai 6 a* 4 d 2 / E_ E + 2 + 8 a x 4 af d* 2 f E_ E+ 2 + 
12ai 6 a* 4 d 2 /*E_E+ 2 + ai 4 aj 6 d* 2 /* E_ E+ 2 + 18 ai 4 af dd* g' E_ E+ 2 + 27 ai 6 af 4 dd* g'* E_ E+ 2 + 9 ai 6 a* 4 d d* / E+ 3 + 
2ai 6 a* 4 dd* /* E+ 3 + 2ai 6 a* 4 d 2 g' E+ 3 + 8 ai 8 af d 2 g'* E+ 3 . 



These polynomials were found using a system of analytical computations Mathematica 4 . 



Appendix 2: 

here we provide the proofs for lemmas stated above. 

1. Due to implicit function theorem, it's sufficient to prove 
that Jacobian of the system © with respect to a set of vari- 
ables (a 3 ,al,X) does not vanish on the straight-line string 
0. Here s > 1. Computing the Jacobian, we have det J = 
2-^(l/| 2 — | 1 2 ) - Substituting the straight-line string solution 
Q to the definitions @, we have: A = |oi|, / = 0, g = {a'j 
at s = 2, at s > 2}. Therefore, det J = {-2|ai| 3 ^ at 
s = 2, at s > 2}. 

2. Let s = 2. In the vicinity of straight-line string we re- 
place g —> a,i + eg', f — > e/, d — » ed, E± — > eE±. Here we 
parameterize the deviation of coefficients from straight-line 
string values by small parameter e and use it as correction 
order counter. Let's prove that 



a2|s ± =o 



t n P n a\ 



^ a* 



(18) 



• Oi • (E+ + /*ci2 + ff'*a 2 + d*a2 2 + ida 2 a 2 ) 



The first two equations in @ with S± = can be rewritten 
as 

02 = 

a 

and it's complex conjugate. Here we have collected all e- 
containing terms at the right hand side. Introducing variable 
e = e/(a* 2 a?) and substituting the expansion 1181 . we obtain 

n>l n>l 

n>l 7i,m>l 
1 j \ * ~n + m + l 7-. 7-1* *2 2 

-fd 2^ 6 p nPmai 0,!. 

n , m > 1 

From here we have the following recurrent relations for P n : 
A = -£+, 

7-1 #■* *2 7-1 * /* 2 j-> 

Pn = -/ Ol P„_i - 3 OiPn-1 

— d a x S ' P n _ m _iP m 
l<m<n-2 

— idaj 2 a? ^ Pn-m-l^mi n > 1. 
l<m<n — 2 

The properties of P„ with respect to i?3 , D4- 
transformations follow from the transformation laws 
1101 for a2, while the other properties of P n can be easily 
proven from their recurrent definition by induction. 

3. Expression 1121 follows immediately from lcg-parame- 
trization 11311 and definition of leg-oscillators given by for- 
mula (5) of Part I. Constraint \3 generates rotation of os- 
cillator variables in plane (a n i,a n 2) and rotation of vec- 
tors e*i,e2 about e-$ in opposite direction, so that the vec- 
tor field a n ± is preserved. xi,2 generate rotations of ez 
and transformations of a n ± by formula 11211 . In this for- 
mula the curve Q as geometrical image is preserved, only 
its parametrization is changed according to il'M . Compo- 
nents of q are proportional to gauge fixing conditions JjjJ: 
q s i = Re a s i = (a s + a_ s + a* + a* a )/4, q S 2 = Re a s2 = 



i(a s + a_ s — a* — a_ s )/4. Evolution, described by Eq- 
transformation: a n ± — > a n ±e~ lnT , is equivalent to a motion 
of marked point O along the curve Q(a), uniform with re- 
spect to parameter 113H : Q(a) — > Q(o + r), or in terms of 

(q,p)- 



COS ST 

— sin st 



sin st 

COS ST 



vectors (q, p) move along ellipses in tangent planes to the 
sphere, see fig|5| In such evolution the points q = move 
along definite curves on the sphere. On such curves the 
vectors q, p are linearly dependent (corresponding ellipses 
degenerate to segments): q x p = 0. The vector q x p ~ 
ia s ± x d* ± is r-independent and always parallel to S3. Thus, 
orbits of q — are defined as zero-level curves of a single 
scalar time-independent function on the sphere: F(Sz) = 
(qxp,e 3 ) = 0. 

4. During the period At = 2-7T all vector fields a n ± are re- 
turned to the initial state. The field a s ± at Ar^ = %2k/s, 
k = 0...S — 1 is returned to the initial state, and at At~ = 
7r(2fc + l)/s, k = 0...S — 1 reverses the sign: a a ± — * —a s ±. 
Singular points q = Re a s ± — are preserved by these trans- 
formations, while the points on the lifting generally are not. 



5. For straight-line string Q(cr) 
ing it to il l 2t we have: 



(cos a, sin a, 0), substitut- 



dcr(— sin a, cos a, 0) exp is(a — V x cos a — V y sin a) 



where V = S3. We see that this expression does not depend 
on V z and gives the same vector fields in the vicinities of 
northern and southern poles. For small V x , y we have in xy- 
plane, omitting unessential constant factors: Six ~ (— i, 1) + 
o(l), (fix ~ (0, l) + o(l) pi± ~ (-1, 0) + o(l) - no singularity; 
a 2 ± ~ (-K - iVy,-iV x + V y ) + o(V x>y ), q 2 ± ~ (-V x , Vy) + 
o(V x , y ), p 2 ± ~ (— Vy, — K) + o(T4,j,) - saddle point; at s > 
2 a s x, q s ±,Ps± ~ (^j;,b) ~ higher order singularity. The 
evolution consists in rotation of vector fields in xy-plane, 
preserving their pole singularities. 

6. The first statement follows from non-degeneracy of sad- 
dle point and general transversality theorems |10| . It also 
can be proven by direct computation: let's consider small 
deviation of solution from straight-line string J3J and as- 
sociated variations of the vector fields at s = 2: a = 
(—V x -iV y , -iV x + V y ) + (Sa x , Sa y ), q = (-14, V y )+(Sq x , Sq y ), 
p = ( — V y , — V x ) + (8p x , Spy). New singular point q = will 
be slightly shifted from the pole: V x — 8q x , V y = —8q y 
(as a result of non-degeneracy of dq x , y /dV x , y ). Consider- 
ing the evolution q — ► gcos2r + psin2r, we have V x = 
Vxo + b cos 4r + d sin 4r , V y = V y o + b sin 4r — d cos At, where 
V x0 = (5q x +5p y )/2, V y0 = (Sq y +5p x )/2, b = (Sq x -5p y )/2, 
d = {Sqy + 6p x )/2, i.e. the saddle point moves along a cir- 
cle near the pole. The same orbit is defined as zero level of 

(qXp, V) = (V X - V x0 f + (Vy - Vy ) 2 -b 2 - d 2 . 

7. From lemma 2 we know, that each monomial in P n con- 
tains at least one operator E(±), additionally n— 1 operators 
from the group (E, d, f,g' and their conjugates) and 2(n — 1) 
operators (ai,af). The numbers of E- and d-operators are 
related: n(E) = n(d) + l. Then we substitute Q to the given 
monomial, expand it to a number of (secondary) monomials 
and normally order each one. Let n(d|), n(/|), n(g'^), n(oi|) 



(2) 

be the numbers of L -raising operators in the secondary 
monomial, n(di), n(fi), n(g[), n(ai|) be the numbers of 

L^-lowering operators. Using the fact, that d-operators 

shift L 2 ' by 4, /, gr'-operators - by 2, by 1 (see Table 1), 

we conclude that in the matrix (: P n :)n 1 n 2 = = JVi | : 

I (2) \ 
Pn '■ Lq 

'i = N2) the given monomial creates non-zero block 
with offsets A T L 2> = 4n(d T ) + 2ro(/ T ) + 2n(g\) + n(ai T ), 



Ai-L, 



(2) 



in(d l ) + 2n(f l ) + 2n{g' l ) + n(a n ) 




N max :>? 




A I 

AiL 



t u 



After summation over all secondary monomials we see 
that the entries in the following triangular block vanish: 
ATI + N2 < N max = A T 4 2) + A^L^ = 4n(d) + 2n(/) + 
2n(g')+n(ai) = 2n(d)+2n(E)-2+2n(/)+2n,(g')+2(n-l) = 
4(n-l). 

8. The implication (i)<S=(ii) is obvious. It's sufficient to prove 
(ii) for one monomial. 

<4 2) < N\o Pln ...o Pt2 o Pll \L^ <N) =Ej Vn -E W3 E JV1 



(4 2) < N\op m \L { 2) = Nn) ■ (4 2) = Nn\ 



(2) 



N2) 



(2) 



(4 2) = 7V2|op l2 |4 2) = JVl) ■ (4 2) = Nl\op tl \L^ < N) 

The summation here is performed over all values of Lq 

(2) 

however, due to the Lq -normal ordering only the values 

(2) 

Lq < N contribute, i.e. the result can be written as a 
product of finite matrices opi n (N)...opi2(N)opn(N). 
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Fig. 7. Spectrum (L Q 2) , A^' ). 



9. The operator a 2 , standing on the right of each : P n : 
in 1160 . annulates all states of Table 2, except of 1 2i2_i^ . 
This state is transferred by af to |2-i). Each monomial of 



: P„ : contains n(E) > 1 operators E± and nj, > L Q 2 ^- 
lowering operators from the group (01, fi, g[, d{), standing 
on the right. Here a± annulates |2-i), other lowering opera- 
tors decrease L 2 ^ by 2,4 and either annulate |2_i^) or trans- 

1 \ (2) 
fer it to 0). The operators E± preserve L and change 

Ag 2 ' by ±1. Figure|7]represents the spectrum of (L Q 2 \ A^). 
Because there are no states immediately above and below 

-1 on this figure, E± 



-(2) 



A 



(2) 



and L. 



annulate |0^ and |2-i). 



(2) 




2, A 



(2) 



Appendix 3: numerical methods. 

The matrix elements of elementary operators (Table 1) in 
\L < N) space are represented by sparse matrices, where 
the typical pattern of non-zero entries is shown on figlHl 
The matrices possess two-level block structure, correspond- 
ing to their (.A3 , Lq )-charge properties. Here large blocks 
have constant A 3 2 ' values, inside them there are smaller 

r (2) 



blocks with constant L . There is an approximate sim- 
ilarity among the blocks. Inside the blocks of the second 
level the matrices possess sparse structure of non-classified 
type ("other!" [H]). 



9 + 



ai 



Fig. 8. Matrix representation of elementary operators 
(TV = 10). 

The most intensive part of the computation is the evalua- 
tion of large polynomials of these matrices. This is done us- 
ing sparse matrix algorithms jll) . where non-zero elements 
of the matrices are stored in condensed form. For sparse 
matrix multiplication two algorithms of these kind are de- 
scribed in .if : pattern multiplication and threshold multi- 
plication. The algorithms compute a product C = A ■ B T , 
where the transposition is computationally expensive oper- 
ation due to associated changes of storage scheme, and is 
performed once for all operators in the Table 1. The pat- 
tern multiplication algorithm is applicable in the case when 
the structure of the product is known in advance and has 
computational cost T\ = O(Ncb), where Nc is a number of 
non-zero entries in C, & is an average number of non-zero 
entries per one row in A, B ("bandwidth"). The threshold 



multiplication algorithm works in the case when the struc- 
ture of the product is unknown and requires T2 = 0(d 2 b) 
operations, where d is the dimension of matrices. Multiplica- 
tion of dense matrices requires Tg = 0(d 3 ) operations, thus 
practically Tj « T2 « T3. We have implemented a com- 
bination of pattern and threshold multiplication algorithms, 
which uses known block structure of matrices and effectively 
takes into account the remaining sparse structure inside the 
blocks. Separability of the computation to the subspaces 
with different values (Az,Q) provides additional optimiza- 
tion. Computing the products prod„ = prod„~i ■ opJ, n s, 
we initialize prodo by the projector to a given (Asi,Qi)- 
subspace, represented as a sparse matrix. Although the in- 
termediate products belong to different (A3, (J)-blocks, this 
initialization automatically cuts off irrelevant contributions, 
considerably reducing the computational time. The local- 
ization of the final result in (A^f, Q f\prod\A3i, QiJ block, 
defined by its charge properties, is useful as consistency test. 
The resulting monomials are injected to a dense matrix ai, 
used to compute a symmetric matrix P 2 by the formula 
1171 . Its eigenvalues are then determined by Householder 
algorithm implemented in the library of dense matrix 
computations NewMat |12| . The computation is iterated for 
different correction orders n < nmax = 6. The scheme of 
main algorithm is given on fig'lHl The necessary input data 
are provided by the following external program modules: 



N,s —A list — H matr \ 



Mathematics 



( NewMat j- 



[for n 


=1; n<=nmax; n++ 


+ 




read_cdeg 




1 


[for 


=0; knmon; i++ j* 



sparse 
dense 



inject 



eigenvalues 



Fig. 9. Scheme of computation. 



12,16,20,24,28 (sector Q = Z^ 2) mod4 = is defined by 
values Nmo&A — 0). There is no A^-dependency in n = 1 
row, because corresponding correction Pi = — E+ commutes 

for this correction arc 



f 2) 

with L , thus non-zero entries of P 



1 f 2) 
located in (A 3 L, 



P 2 A^'L, 



(2) r (2)\ 



/ 



diagonal blocks. When N 



is increasing, new blocks appear, which give contributions to 
higher regions of P 2 , leaving intact already filled eigenval- 
ues. Starting from n = 2, L -nondiagonal contributions 
appear, leading for some eigenvalues to iV-dependence sat- 
urated at large N. Other eigenvalues, as explained earlier, 
correspond to the states annulated by 02 and remain to be 
constant, integer and degenerate. Spectra in other Q-sectors 
show similar behavior. 
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Fig. 10. P 2 /2 it as a function of n, ^-parameters. 



Time and memory requirements of the algorithm are pre- 
sented by the following statistics: 



N 


12 


16 


20 


24 


28 


subdim 


69 


258 


890 


2851 


8567 


req. memory 


1.4Mb 


4.4Mb 


23Mb 


147Mb 


1.1Gb 


comp.timc 


2 sec 


44 sec 


14 min 


4 hours 


66 hours 



Here subdim is the dimension of largest (A3 = 0, Q = 0) 
block in dense matrices, which at large N mainly defines the 
amount of required memory. The computation is performed 
on HP 2GHz Linux PC. 



list - computes for the given N, s the spectrum of 
(I/q S ' , A^) (see fig|7J and corresponding basis in the 



space 



-00 
^0 



<N); 



• matr - computes the matrix elements of the elementary 
operators in this basis; 



Mathematica script, which substitutes the defi- 
nitions of elementary operators in the polynomial ex- 
pressions of Appendix 1, expands them to normally or- 
dered monomials and represents the result in the form 
of coefficient-degree lists. The number of monomials for 
each order n = 1..6 is given in the following table: 



n 


1 2 3 4 5 6 


num. of monomials 


1 4 19 82 297 922 



The determined spectrum of mass (P 2 /2tt, S = 0, Q = 0) 
as a function of the correction order n and cutoff param- 
eter in Lq 2 ' < iV is presented on fig llOl The spectrum 
is subdivided to rows with fixed n. The points in each 
row show, in the order up-to-down, the spectra for N = 



